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Abstract 

We estimate analytically the second-order cosmic microwave background temperature anisotropies at 
the recombination epoch in the squeezed limit and we deduce the contamination to the primordial local 
non-Gaussianity. We find that the level of contamination corresponds to = 0(1) which is below the 
sensitivity of present experiments and smaller than the value 0(5) recently claimed in the literature. 
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1 Introduction 

Primordial non-Gaussianity (NG) of the cosmological perturbations has become a crucial aspect of both 
observational predictions of (inflationary) early universe models and of present and future observational 
probes of the Cosmic Microwave Background (CMB) anisotropies and Large-Scale Structures (LSS) of 
the universe. The main motivation is that detecting, or simply constraining, deviations from a Gaussian 
distribution of the primordial fluctuations allows to discriminate among different scenarios for the gener- 
ation of the primordial perturbations. A non-vanishing primordial NG encodes a wealth of information 
allowing to break the degeneracy between models that, at the level of the power spectra, might result 
to be undistinguishable [1]. It is impressive that a future detection of a high level of primordial NG 
could rule out the standard single-field models of slow-roll inflation, which are characterized by weak 
gravitational interactions and thus predict tiny deviations from Gaussianity [21 [3j. Statistics like the 
bispectrum (three-point correlation function) of CMB anisotropies [1], and higher-order statistics, and 
various observational probes of the LSS [5], can be used to assess the level of primordial NG on various 
cosmological scales. 

It has now been fully appreciated the importance of distinguishing different shapes of the primordial 
bispectrum, i.e. the dependence of the bispectrum on different triangular configurations formed by the 
three wave- vectors over which, in Fourier space, the primordial fluctuations are correlated [6]. Different 
(classes of) inflationary models give rise to unique signals with specific triangular shapes, which thus probe 
different aspects of the physics of the primordial universe, in particular they can probe different kinds 
of interactions in the inflaton sector which are at the origin of the primordial non-Gaussian signal. For 
example, models in which non-linearities develop outside the horizon during inflation, or immediately after 
inflation, generate a level of NG which is local, as the NG part of the primordial curvature perturbation 
is a local function of the Gaussian part, being generated on superhorizon scales. This is generally the 
case when a light scalar field, different from the inflaton, contributes to the curvature perturbation (for 
example, multi-field models of inflation [7j, curvaton [HI [9] and inhomogenoeus reheating models |10 | lll [ 
112]). In momentum space, the three point function, or bispectrum, arising from the local NG is dominated 
by the so-called squeezed configuration, where one of the momenta is much smaller than the other two 
(k\ <C &2 ~ ^3) an d it is parametrized by the non-linearity parameter "Equilateral" NG, peaking 

for k± ~ &2 ~ &3 is generally predicted by inflationary models where the inflaton field is characterized 
by derivative interactions so that the correlation between the modes peaks when all the modes cross the 
horizon during inflation at the same time (see, e.g., inflationary models with non-standard kinetic term 
|13j . or DBI models of inflation j!4j). For other interesting NG shapes/models (such as "orthogonal" 
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[15] and "folded" (or flattened) shapes JT31 HH [T71 [TBI IS] ) we refer the reader to the review [20] . Present 
limits on NG are summarized by the WMAP measurements -10 < < 74 and -214 < f^ 1 < 266 at 
95% CL [21] (see [22] for CMB constraints on other types of NG). The Planck satellite [23] is expected 
to improve the sensitivity to the non-linearity parameters by an order of magnitude. In particular, for 
the local case, it is expected to be sensitive to a value of the nondinearity parameter as low as /^l = 3 at 
la (accounting also for polarization information; for the equilateral case /^ L uil = 30) [Ml [251 [261 [23 [28] . 

However there are many sources of NG in the CMB anisotropics beyond the primordial one, including 
systematic effects and astrophysical contamination. Given the deep implications that a detection of 
primordial non-Gaussianity would have, it is crucial that all possible sources of contamination for the 
primordial signal are well understood and under control. In fact any nondinearities can make initially 
Gaussian perturbations non-Gaussian. Such non-primordial effects can thus complicate the extraction 
of the primordial non-Gaussianity. One of their main effects could be that of "mimicking" a three-point 
correlation function that is similar in shape to the primordial one, thus inducing a bias to the estimator 
of primordial NG used. Therefore we have to be sure we are not ascribing a primordial origin to a signal 
that is extracted from the CMB (or LSS) data using estimators of non-Gaussianity when that signal 
has a different origin. Moreover we must always specify which primordial non-Gaussianity we study the 
contamination from the non-primordial sources (e.g. primordial non-Gaussianity of "local", "equilateral" 
or "folded" shape). In that respect the reasonable question to ask is which kind of primordial non- 
Gaussianity a given secondary effect can contaminate most, wether it is of the so-called local type or of 
the equilateral type. 

As far as CMB anisotropics are concerned, well-known examples of contamination are the NG pro- 
duced by (non-linear) secondary anisotropics, like, e.g., the Sunyaev-Zeldovich effect, gravitational lens- 
ing, Integrated-Sachs- Wolfe (or Rees-Sciama) (see, e.g. [29] for a review on CMB secondary anisotropics). 
Up to now it has been found that the most serious contamination to the local type of primordial NG is the 
the coupling between the Integrated Sachs- Wolfe/Rees-Sciama effect and the weak gravitational lensing 
[30] . In the detailed analysis of [HI [32], [33] it is shown that the ISW/Rees-Sciama-lensing bispectrum 
produces a bias of /nl — 10 to local primordial non-Gaussianity, while the bias to equilateral primordial 
non-Gaussianity is negligible (see also [MIES])- Other secondary effects have a little impact on the (local) 
primordial NG with |/nl| < 0.7 (for a summary of the level of contamination see the reviews [361 14"] and 
references therein). 

Another relevant source of NG comes from the physics operating at the recombination epoch. The 
interest in these sources of contamination is due to the fact that they actually bring previously unknown 
effects coming from the non-linearities in the Boltzmann equations. A naive estimate would tell that 
these non-linearities are tiny being suppressed by an extra power of the gravitational potential. However, 
the dynamics at recombination is quite involved because all the non-linearities in the evolution of the 
baryon-photon fluid at recombination and the ones coming from general relativity should be accounted 
for. This complicated dynamics might lead to unexpected suppressions or enhancements of the NG at 
recombination. A step towards the evaluation of the three-point correlation function has been taken in 
Ref. [37] where some effects were taken into account in the squeezed triangle limit, corresponding to the 
case when one wavenumber is much smaller than the other two and was outside the horizon at recom- 
bination. Refs. |38[ [39"1 |4"0"1 14T1 142] . present the computation of the full system of Boltzmann equations, 
describing the evolution of the photon, baryon and Cold Dark Matter (CDM) fluids, up to second order in 
the perturbations. These equations allow to follow the time evolution of the CMB anisotropics at second 
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order on all angular scales from the early epochs, when the cosmological perturbations were generated, 
to the present time, through the recombination era. These calculations set the stage for the computation 
of the full second-order radiation transfer function at all scales and for a a generic set of initial conditions 
specifying the level of primordial non-Gaussianity Q. 

There has been some debate about the level of contamination that the non-linear evolution at re- 
combination can generate. A first work, partially implementing numerically the second-order Boltzmann 
equations in the small scale limit, claimed a very high level of contamination to the local-type /nl [43] . 
However it has been shown in Ref. |44| that the source of this effect, the small-scale non- linear evo- 
lution of the second-order gravitational potential at recombination, actually contaminates mainly the 
equilateral-type NG. A contamination of /^^(equil) ~ 5 for an experiment like Planck was found [44] . 
which is actually low w.r.t. the expected Planck sensitivity to this kind of NG. The contamination to 
the local NG from the small scale evolution of the second-order gravitational potential is also very low, 
/nl^O 00 ) — 0-3 [S] for an experiment like Planck. In Ref. [45] another numerical implementation of 
the second-order Boltzmann equations has been performed to compute the contamination from specific 
terms appearing in the second-order radiation transfer function. These are the terms that can be written 
(in the Poisson gauge adopted in Ref. [45]) as products of the first-order perturbations (in the form, e.g., 
(c^ 1 )) 2 , where $W is the linear gravitational potential). The contamination to local NG was found to 
be negligible for an experiment like Planck being (/^^(loc)! < 1 [45] . 

Recently, Ref. [46] analysed numerically the second-order evolution of the photon-baryon system at 
recombination claiming an overall contamination /^l^O 00 ) = 5 up to a m axi mum multipolc f m ax — 

2000. 

Such a contamination must therefore come from terms in the second-order radiaton transfer function that 
are "intrinsically" of second order. The authors of Ref. [46] suggest that the main candidate is the second- 
order monopole perturbation (6sw = ^oo /^ + ^^ 2 ^)> which is the usual term appearing in the CMB 

(2) 

anisotropies due to the intrinsic photon energy density fluctuations Aq and the gravitational redshift 
due to the potential, both evaluated at recombination. In fact such a term on large scales reduces to 
the Sachs- Wolfe effect. There has been some debate on how such a second-order contribution is able to 
produce a contamination /^^(loc) = 5 (see, e.g., [471 08]). For example it is not clear what is the physical 
mechanism generating such a squeezed signal: if it occurs due to an integration over the terms (squared in 
first-order perturbations) that source the "intrinsically" second-order monopole, then, naively, one does 
not expect that a squeezed signal can be generated, because such a contribution would be associated to 
a causal evolution on small scales. 

The goal of this paper is to present a transparent computation of the bispectrum in the squeezed limit 
through a convenient coordinate rescaling \37 \ I49 [ [50]. To understand such a rescaling, it is important to 
recall what is generally the origin of a squeezed non-Gaussian signal: typically the local-form bispectrum 
is generated when short-wavelength fluctuations are modulated by long-wavelength fluctuations. In 
particular we will focus on the temperature anisotropies at recombination when the long wavelength 
mode is outside the horizon, but observable at the present epoch. Thus, the effect of the long wavelength 
mode imprinted at recombination can be described simply by a coordinate transformation. In this way 
we can describe in a simple way the coupling of small scales to large scales that can generally produce the 

1 Of course, for specific effects on small angular scales like Sunyaev-Zeldovich, gravitational lensing, etc., fully 
nonlinear calculations would provide a more accurate estimate of the resulting CMB anisotropy, however, as long 
as the leading contribution to second-order statistics like the bispectrum is concerned, second-order perturbation 
theory suffices. 
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local form bispectrum. A similar cross-talk between large and small scales gives rise to the ISW-lensing 
cross-correlation bispectrum. 

The plan of the paper is the following: In Sec. II we show that a simple coordinate rescaling reproduces 
correctly the second-order expressions for the temperature anisotropies in the squeezed limit. In Sec. Ill 
we take advantage of this coordinate rescaling to estimate the bispectrum from the CMB anisotropies 
at recombination in the squeezed limit and the corresponding contamination to the local primordial 
non-Gaussianity. Sec. IV contains our conclusions. 



2 Second-order CMB anisotropy at recombination 

To motivate the fact that the rescaling of coordinates (that we will soon introduce) correctly captures 
the physics of the CMB anisotropies at second order in the squeezed limit, we will first show that 
such a rescaling applied to the linear equations properly reproduces the second-order sources found in 
Refs. \36\ [38l I39| . We first recall the Boltzmann equations at first-order in perturbation theory for the 
CMB anisotropies. Since this subject is rather standard, we refer the reader to standard books for more 
details [51] (see also Refs. [361 EH]). O ur starting metric is 



ds 2 = a 2 (7]) 



2<I> j 2 i -2* i 2 

-e ar] + e dx 



(1) 



where a{rj) is the scale factor as a function of the conformal time rj, and we have neglected vector and 
tensor perturbations. The equations of motion of the first two moments of the Boltzmann equations for 
the CMB photons are 

A 

0, (2) 



(3) 



where IP- 7 is the photon quadrupole moment and r' = —n e aj-a is the differential optical depth in terms 
of the average number of electron number density n e and Thomson cross section ax- Here the primes 
indicate differentiation with respect to rj and di differentiation w.r.t. x l . The two equations above are 
complemented by the momentum continuity equation for baryons, which can be conveniently written as 

R 



4 1 * 



(4) 



where we have introduced the baryon-photon ratio R = 3pf,/(4p 7 ) and we have indicated by H = a' /a 
the Hubble rate in conformal time. 

Eq. (|4|) is in a form ready for a consistent expansion in the small quantity t~ 1 which can be performed 



in the tight-coupling limit. By first taking v^ 1 
left-hand side of Eq. one obtains 



at zero order and then using this relation in the 



Such an expression for the difference of velocities can be used in Eq. 
for the photon velocity in the limit of tight coupling 



,(1)* 



R 



1 A (l),i 
l + R^ Al + R 



0. 



(5) 

to give the evolution equation 

(6) 
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Notice that in Eq. (|6|) we are neglecting the quadrupole of the photon distribution n^ 1 ^- 7 (and all the 
higher moments) since it is well known that at linear order such moment(s) are suppressed in the tight- 
coupling limit by (successive powers of) 1/r' with respect to the first two moments, the photon energy 
density and velocity. Eqs. ([2]) and (|6j) are the master equations which govern the photon-baryon fluid 
acoustic oscillations before the epoch of recombination when photons and baryons are tightly coupled by 
Compton scattering. 

In the tight-coupling limit the contribution to the CMB anisotropics from recombination is given by 

e (1) (k,n,i,„) = Qa^ + ^W+vW -n) . (7) 

We are now interested in computing the second-order CMB anisotropies in the squeezed limit, i.e. when 
one of the momenta is much smaller than the other two, k\ <C ki ~ k%. In particular, we will be interested 
in the case in which the smallest momentum is at most equal to the wavenumber entering the horizon at 
equivalence, k\ < k eq . 

Instead of solving the complicated network of second-order Boltzmann equations for the CMB temper- 
ature anisotropies, we use the following trick. As the wavenumber k\ <■ k eq corresponds to a perturbation 
which is almost larger than the horizon at recombination and the evolution in time of the correspond- 
ing gravitational potential is very moderate (one can easily check, for instance, that changes 
its magnitude by at most 10% during the radiation epoch for k\ = k cq ), we can absorb the large-scale 
perturbation with wavelength ~ k^ 1 in the metric by redefining the time and the space coordinates as 
follows. Let us indicate with <Iv an d the parts of the gravitational potentials that receive contributions 
only from the large-scale modes k\ < A; cq (see, e.g., (49] ) . If the scale factor is a power law a{r\) oc r] a 
(a = 1 and a = 2 for the period of radiation and matter domination, respectively), we can perform the 
redefinitions 

a^e^dr/ 2 = r^V^dr/ 2 = rf a drf = a 2 (r?)drj 2 =j> 77 = e Tf=% , (8) 

and 

a 2 (r ? )e- 2 *Mx 2 = T^e^dx 2 = rj 2a e~ 2 ^ dx 2 = a 2 (rj)dx 2 => x = e *T=*« e -*'x . (9) 

In particular, the combination 

kjj = e *«+*t kr) , (10) 

where k and k are the wavenumbers in the two coordinate systems. Obviously, if one wishes to account 
for the fact that at recombination the universe is not fully matter-dominated, one should perform a 
more involved coordinate transformation which will eventually depend also on the parameter R. To test 
the goodness of this procedure, let us consider for instance the equation of motion of the gravitational 
potential x I / k(??) at second-order during the radiation epoch. From the full second-order equations (B.13) 
and (B.14) of Ref. [36], we see that this equation reads in the squeezed limit 

< 2) " + m^' + clkH? = _ 4c? + 41)) t 

where in the source in the r.h.s. we have implicitly assumed an integration over the momenta ki and 
k2 with the corresponding Dirac delta to ensure momentum conservation. We will adopt this convention 
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from now on. In fact, for a bispectrum computation, it is the kernel of the source term that matters, 
and from this point of view we can identify the wavenumber ki as a long wavelength modulation of the 
perturbation in the momentum k2, so that k = ki + k2 — k2. 

It is reassuring that this equation can be simply obtained by rescaling the coordinates, following (JSj) 
and ([9]) . Just consider the (linear) equation of motion [511 [39] 

$(i)" + 4 W (i)'_ c 2 V 2$(i) = o, (12) 

whose solution is 

, T ,(i)^ - o lTr (i)/ m (-c B fcr/cos(c 8 fcTy) +8in(c 8 fcTy)) 
k K ,) k y ' (c s k7])3 ' K } 

Applying © and © to Eq. flT2|) we obtain 

*" + m$' - c 2 e 2(<Iv+ * £) V 2 1- = . (14) 

Now, expanding the perturbations into a first- and second-order parts, VP = + $( 2 ) /2, and identifying 
in Fourier space the long wavelength mode with ki, one can easily check that Eq. (|14h exactly correponds 
to Eq. (HH). The solution of Eq. ([II]) is 



* (2) (j7) = 3'I^ 2) (0) (-^ fc? ? cos (^ fc? ?) + sax{c a krj)) 



[c s kr]) 3 

+ g 3c s kricos(c s kri) + (-3 + (c s k V ) 2 ) sm{c s kri) _ 

(c s kr)) 3 kl k2 

Again, to have a further check of our procedure, we point out that this solution can be found by employing 
the coordinate rescalings ([8]) and ([9]). Taking the expression fjl3f) and assuming that ^ ~ we find 



*kfa) = 3%(0) 



-c s k 2 rjcos(c s k 2 r]) + sm(c s k 2 r]) 
(c s k 2 r]) 3 



-e k i c s /c2?? cos I e k i c^f? ) + sin I e k i c s k 2 rj 
3%(0p / ~ J — ~ —■ (16) 



2* 



(i) 



c s k 2 r] 



Expanding at second-order one, one recovers exactly the expression (jlSjl 2 !. 

We have also verified that using this procedure applied to Eqs. ([2]) and ([6]) one reobtains the cor- 
responding second-order Boltzmann equations already computed in Eqs. (205) and (206) of Ref. |36| 
expressed in terms of the bolometric temperature, i.e., in terms of the variable [A^ — (A^) 2 ]. The bolo- 
metric temperature is defined as (see, e.g., [H] 05]) (T/T) 4 = A, where A(x, n, r) = J dpp 3 f/(f dpp 3 f) 



Notice that here we are using the short-hand notation kr\ = e k i k i kr\ which simplifies the computation 
giving nonetheless the correct result. Notice that the coordinate rescaling also transforms the gravitational potential 
v E'k(0). This change however is reabsorbed in the physical initial conditions. 
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is the (normalized) brightness function, / being the photon distribution function. By Taylor expanding 
this definition one finds 

A (2) = 4@ (2) + 16 ( (1)) 2 = 40 (2) + ^ A (l)^ 2 j (17) 

where, in the notations of Refs. [Ml [36] . 

t = Te e = t , ^ = ew + ie( 2 ' + i(e( 1 )) 2 , (is) 

is the bolometric temperature so that one recovers the temperature variables used in [45, 36J (notice that, 
despite using the same symbol G this temperature is different from the brightness temperature defined, 
e.g., in Eq. (4.1) of [46]) . Since the existing literature usually computed the CMB bispectrum in terms of 
the bolometric temperature (e.g., [451 146]). our results for the CMB bispectrum will refer to this variable. 
Let us show in some details how the coordinate rescaling method works for Eqs. ([2]) and ([6]). Under the 
coordinate rescalings dH) and ([9]) Eq. ([2]) transforms as 

4e , 00 + ^+*^4-4* / = 0. (19) 

Expanding this equation at second-order in the perturbations (and expressing the photon velocity in 
terms of Aoo and \I/) one finds 

49$' + \divf* - W' {2) = -2(d>« + - 46$') . (20) 

This equation can be confronted with Eq. (205) of Ref. [36] 

AS' + - W {2Y = S A , (21) 

where the source term is given by 

S A = (a$ 2 )'-2 ( <I>M + ^ 

16 8 / H m, I^aF' 



(22) 



It is easy to check that in the squeezed limit, ki <C &2 ~ Eq. (|2T|) exactly coincides to Eq. ([20]) when 
written in terms of the bolometric monopole [A^ — (A^) 2 ]. Similarly one can check that the squeezed 
limit of Boltzmann equation for the second-order velocity equation of the photon-baryon system can be 
obtained in the same way. Let us start from the (linear) Eq. ([6|). By applying the coordinate rescalings 
([E]) and © we obtain 

< + « YTR^ + \ e ^ ] 1+1 + e( * + * ) = ' (23) 

which, expanded at second-order in the perturbations, becomes 

v ? 1 ' + n r^ v ? 1 + tttSt + * (2),i = -^rrrm^ {1) + * (1) ) 4 ©oo - 2(* (1) + tfW)** 1 ^ . (24) 

1 + it 4 1 + Jrt 2(1 + tx) 



This equation corresponds to the squeezed limit of the velocity continuity equation computed directly at 
second order in Ref. [36] (see Eq. (209)) 

ft 1 a( 2 ^ 

v^> + H^vW + \=^ + ^ = 4 , (25) 



where 



5v " "iTfl) 9j ^ - 2w *- 2W iTb w + 2 (TT1?F 00 7 

+ 4(1 + J R)2 V A oo J + 3(l + i?)^ °^ + l + R 

_ 2(c&( 1 ) + - 1 ($(D + V^A^ 

2(1 + i?) 

^ _«<„(i)2 _ 3 g A (1) f - - I^L^ C2fil 



1 + i2 7 21 + i? u Ml + J R 7 Al + R 



Here u/j is the pure second-order metric perturbation of the (0 — i) metric tensor (which therefore does 

(2) 

not appear in the first-order equation) while H\j is the second-order quadrupole moment of the photons. 
If one takes the source term S v , Eq. (126f) . in the squeezed limit and for R = 0, one recovers exactly 
Eq. ([HP once Eq. (f25j) is expressed in terms of the bolometric temperature [Aqq — (A^) 2 ]. Notice 
that the second-order velocity continuity equation of the photon-baryon system in the squeezed limit is 
recovered only when k\Rrj rcc < 1, i.e. when the timescale of the modulating long wavelength mode is 
much bigger than the typical timescale of the collision term in the tight coupling limit (in particular for 
]?< 1). This does not come as a surprise, the coordinate transformation we are adopting is not fully 
exact when R is not much smaller than unity. 



3 The bispectrum at recombination in the squeezed limit 
and the contamination to the primordial local NG 

Since we will be concerned with a signal-to-noise ratio (S/N) dominated by the maximum multipole a 
given experiment can reach, ^ max S> 1, we can use the flat-sky approximation [52] and write for the 
bispectrum 

(a(f 1 )a(f 2 )a(f 3 )) = (2tt) 2 5^(£ 1 + £ 2 + £ 3 )B(l 1 ,£ 2i e 3 ) . (27) 

Inspecting Eq. f)18|) . the bispectrum gets two contributions. One from the intrinsically second-order term 
6( 2 ) and one from the term (0^)) 2 /2. 

The first contribution can be calculated through the rescaling of the coordinates. We make the 
simplifying assumption that rj rcc S> r] eq in such a way that the coordinate transformations ([8]) and ([9]) 
can be performed in a matter-dominated period, that is we take a = 2 and 

x ^ e -5«v/3 X) k _> e 5*</3 k) Viec ^e**/ 3 r} iec , (28) 
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is the transformation for modes which were outside the horizon at recombination, but are subhorizon at 
the time of observation. It is therefore enough to consider the rescaling acting on 

a[i) = / 1>7 ** (0) At(£ ' fc2) ' (29) 

where At(£, k z ) is the radiation transfer function and D iec = (770 — rj rcc ) is the distance to the surface of 
last scattering. We consider only the contribution at recombination and therefore 



MA k z ) = *S{J(P) 2 + £ 2 / D ?cc rirec) , S = Qa$ + *W + vW • n) /<&«(0). (30) 



r)2 W v y 7 rec ' /rec; ' \ 4 

-^rec V^ 1 / r]=r] Iec 

Here fc 2 is the momentum component perpendicular to the plane orthogonal to the line-of-sight and the 
superscript reminds us that in the expression for k we have to set P = £/D*, where P is the component 
perpendicular to the line-of-sight. As we have learnt in the previous section, the rescaling (|28|) reproduces 
the second-order sources and the corresponding solutions starting from the first-order expression of the 
temperature anisotropy. Therefore, we have to think of the expression (130)) as a function of the rescaled 
coordinates in such a way that an expansion of the long mode will give the exact second-order quantity 
in the squeezed limit. Notice that the rescaling (|28j) changes also the gravitational potential 



41) _> e -», e/3k , (31) 

and therefore a(£) is not subject to any rescaling if one takes the Sachs- Wolfe large-scale limit in which 
the source S reduces to 1/3. 

Following Ref. |44j . we mimic the effects of the transfer function on small scales as 

At (£, k z ) = 4- VV^I/M" (32) 

D 

i.e. a simple exponential and a normalization coefficent a to be determined to match the amplitude of 
the angular power spectrum at the characteristic scale I ~ I* = k*D rec . As shown in Ref. the values 
— 750 and a ~ 3 are able to account for the combined effects of "radiation driving", which occours 
at t > £ oq ~ 160 and boosts the angular power spectrum with respect to the Sachs- Wolfe plateau, and 
the effects of Silk damping which tend to suppress the CMB anisotropics for scales t > £^ ~ 1300. The 
combination of these effects produces a decrease in the angular power spectrum from a scale I* ~ 750. 
The choice of the exponent 1.2 derives from the study of the diffusion damping envelope in Ref. )53j . 
The power spectrum in the flat-sky approximation is given by 

Hh)a(T 2 )) = (27r) 2 5^(h + h)C{l x ) , (33) 

with 

C(£) = fe/ dk z \A T (£,k z )\ 2 P(k). (34) 
The exponential of the transfer function allows to cut off the integral for k — k„ and one finds [44] 

C(£)ca^ £ f 3 e-W^\ (35) 
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which holds for £ 3> £* and we have used the amplitude of the primordial gravitational potential power 
spectrum computed at first-order 

($(i)(k 1 )$( 1 )(k 2 )) = (2^) 3 <5( 3 )(k! + k 2 )P(fci) , P(k) = A/k 3 (36) 

with amplitude A = 17.46 x 10~ 9 . To compute the bispectrum, we go to the squeezed limit £\ <C £ 2 , £3 (or 
k\ <S k 2 , k 3 ). In this case ^ acts background for the other two modes. One can therefore compute 
the three-point function in a two-step process: first compute the two-point function in the background of 
an d then the result from the correlation induced by the background field. Following Ref. |37j (see 
also [54J) and using the Sachs- Wolfe limit for the multipole £\, this procedure leads to 

{a(£ 2 )a(£ 3 )) w = {a(£ 2 )a(£ 3 )) + 5a(£ 2 + £ 3 )C{£ 2 ) d jf . (37) 

*kj d In £ 2 

The bispectrum from the intrinsically second-order term therefore reads 



B B m{h,t 2 ih) = (a(£i)(a(£ 2 )a(£ 3 ))) = (2vr) 2 6®(£i + £ 2 + £3) 5 C{h)C{l 2 ) ^M^ZM . (38) 

Notice also in all these computations one can safely neglect the rescaling of the conformal time at 
recombination r/ rec up to terms r] Tec /r]o <C 1. 

Next we compute the contribution from the term (Q^) 2 /2. In the squeezed limit it reads (see also 
Refs. 05l[36]) 

5i (e{1))3 (4,4,4) = (2vr) 2 6^(h + £ 2 + £ 3 ) 2 C{£ X )C{1 2 ) . (39) 
The total bispectrum therefore is 

B iec (h,£2,£3) = (2n) 2 5^(£ 1 + £ 2 +£ 3 )C(£ l )C(£ 2 ) 

Notice that if we consider the second-order CMB anisotropics on large-scales with all the modes on super- 
Hubble scales at recombination (i.e. with all wavenumbers (kiT] Tec C 1), but inside the horizon today), 
then in the Sachs- Wolfe limit H 

f = ^ + ^(*S) 2 . (41) 

The term ($rec 2 /18) corresponds in this limit exactly to the (G^ 1 ) 2 /2) term in Eq. (|18|) . The result (|4ip was 
first obtained in Ref. [55] (see also the discussion that followed in Ref. [50J ) . This additional contribution 
to primordial non-Gaussianity from the second-order gravitational potential in ($ rec = ^rcc + <J?rec/2) 
and the term (3>rec (see Refs. [55^ H9]) can contribute to a contamination of the primordial local 
non-Gaussianities only on the very largest scales, 2 < £ < 200. Therefore they actually lead to a small 
(negligible) contribution to the total level of contamination (see the quantity defined below). 

3 This expression is valid only when all the wavelength modes are outside the horizon at recombination, and, 
e.g., when written as a convolution in Fourier space, it does not hold for an external wavenumber kii] Tec <C 1 with 
the two internal wavenumbers kj, and kg inside the horizon. 



2 + 5 



din [£ 2 2 C(£ 2 )] 
din £ 2 



(40) 
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Our goal now is to estimate the level of degradation that the NG from recombination in the squeezed 
limit causes on the possible measurement of the local primordial bispectrum. A rigorous procedure is to 
define the Fisher matrix (in flat-sky approximation) as 

p.. — f Si, d 2 lod 2 ^ <5 (2) it + to + &) ^(^2,13)^(11,12,13) (4 ) 

where % (or j)= (rec,loc), and to define the signal-to-noise ratio for a component i, (S/N)i = 1/J F^ 1 , 
and the degradation parameter = FuF^ 1 , due to the correlation bewteen the different components 

measure the contamination to the primordial bispectra one can define 
that effective non-linearity parameter which minimizes the x 2 defined as 



X 2 = / d 2 hd 2 £ 2 d 2 £ 3 6® (£ + £ 2 + £3) 



B ioc(h, £2, £3; f^l = 1) - B rec (£i,£ 2 ,£3) 



SC{h)C{l 2 )C{h) 



to find 



con 



rcc,loc 



(43) 



(44) 



r 



-^loCjloc 

In multipole space the bispectrum induced by a local primordial NG in the squeezed limit is given by 

(45) 



Bx oc {h,h,h) = -Qfh°l [C{£i)C{£ 2 ) + cycl 
Performing the integrals in the multipoles and using Eq. ()35|) . we find 



1 5 

/•con ^ 1 

/nl - "g + y 2 



1 + 0.75 



1.2 



1 + 0.75 



1.2^ 



~ 0.91, 



(46) 



where we have taken £ max — 

2000, 4 = 750 and 

£mm — 1200, see Ref. 

Let us also remark that in Refs. [44 \ l3"o] the non-linear evolution on small scales (for all the wavenum- 
bers (ki7] rcc ^> 1) at recombination) of the second-order gravitational potential also produces a small 
contamination to the local non-Gaussianity of the order of = 0.3. This contribution should be added 
to Eq. ([MD. 



4 Conclusions and remarks 

In this paper we have analytically estimated the level of non-Gaussianity produced by the non-linear 
evolution of the photon-baryon system at recombination in the squeezed limit. While a contamination 
0(5) was numerically obtained in |46j . we do find that the total contamination to the primordial local 
non-Gaussianity is smaller and not within the reach of present experiments. Our main goal was to provide 
a clear and simple way to understand the physical origin of such a contamination. We have reached the 
goal using a simple rescaling of the local coordinates of a perturbed FRW universe, which explains in a 
very transparent way the large-scale modulation of the perturbations which is generally at the origin of 
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a squeezed non-Gaussian signal. Let us make some brief remarks about the validity of our results. It is 
clear that they are valid up to terms V3>j^ ), which vanish in the exact squeezed limit, k\ — > 0. 

Of course such contributions will arise for the perturbation modes ki that were not much outside the 
Hubble radius at recombination, e.g. for k\ ~ fc eq . For example the initial condition for the second-order 
gravitational potential would get a small correction (see Eq. (269) of |36| ) from the initial condition (with 
vanishing primordial non-Gaussianity) 

4 2) (0) = H F(kl ^ 2 ' k) < ) (0)< ) (0), (47) 

where the function i ? (ki,k2,k) would vanish in the exact squeezed limit (see its expression in Eq. (242) 
of [36]. It is easy to check that, if for example, one takes k\ = k cq and &2 = 15fc e q (the latter corresponding 
to £2 = 2000) one finds that the contamination to the primordial local NG induced from the correction of 
the F term in Eq. (|47p is [/^(loc)! ~ 0.1. Indeed this must also be intended as an upper limit, since this 
correction does not fully correlate with the local type /^l- This example is to show that the corrections 
coming from the perturbation modes k\ = k eq is actually subdominant. 
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Note added 

When completing this work, we have become aware of a similar work by P. Creminelli, C. Pitrou and F. 
Vernizzi. Our results, when overlap is possible, agree with theirs. We thank them for useful correspon- 
dence. 

References 

[1] N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Phys. Rept. 402, 103 (2004) 
[2] V. Acquaviva, N. Bartolo, S. Matarrese and A. Riotto, Nucl. Phys. B 667, 119 (2003) 
[3] J. Maldacena, JHEP 0305, 013 (2003). 

[4] For a recent review, see, e.g., M. Liguori, E. Sefusatti, J. R. Fergusson and E. P. S. Shellard, Adv. 
Astron. 2010, 980523 (2010) 

[5] For example, see [3] and C. Carbone, L. Verde and S. Matarrese, Astrophys. J. 684, LI (2008) 

[6] D. Babich, P. Creminelli and M. Zaldarriaga, JCAP 0408, 009 (2004) [arXiv:astro-ph/0405356l . 

[7] N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 65, 103505 (2002); F. Bernardeau and 
J. P. Uzan, Phys. Rev. D 66, 103506 (2002); F. Vernizzi and D. Wands, JCAP 0605, 019 (2006). 



13 



D. H. Lyth, C. Ungarelli and D. Wands, Phys. Rev. D 67 (2003) 023503. 
N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 69, 043503 (2004). 
T. Hamazaki and H. Kodama, Prog. Theor. Phys. 96 (1996) 1123. 

G. Dvali, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 69, 023505 (2004); L. Kofman, 
|arXiv:astro-ph/0303614l 

G. Dvali, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 69 (2004) 083505. 

X. Chen, M. x. Huang, S. Kachru and G. Shiu, JCAP 0701, 002 (2007) [arXiv:hep-th/0605045l . 

M. Alishahiha, E. Silverstein and D. Tong, Phys. Rev. D 70, 123505 (2004) [arXiv:hep-th/0404084] . 

L. Senatore, K. M. Smith and M. Zaldarriaga, JCAP 1001, 028 (2010) [arXiv: 0905 .37461 [astro- 
ph.CO]]. 

R. Holman and A. J. Tolley, JCAP 0805, 001 (2008) |arXiv:0710. 13021 [hep-th]]. 

N. Bartolo, M. Fasiello, S. Matarrese and A. Riotto, JCAP 1008, 008 (2010) [arXiv: 1004.08931 
[astro-ph.CO]]. 

P. D. Meerburg, J. P. van der Schaar and P. S. Corasaniti, JCAP 0905, 018 (2009) [arXiv:0901.4044l 
[hep-th]]. 

P. Creminelli, G. D'Amico, M. Musso, J. Norena and E. Trincherini, JCAP 1102, 006 (2011) 
[arXiv: 101 1.30041 [hep-th]]. 

X. Chen, Adv. Astron. 2010, 638979 (2010) }arXiv:1002.14l6l [astro-ph.CO]]. 

E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) [arXiv: 1001 .45381 
[astro-ph.CO]]. 

J. R. Fergusson, M. Liguori and E. P. S. Shellard. larXiv: 1006. 16421 [astro-ph.CO]. 



See http://planck.esa.int/ 



E. Komatsu and D. N. Spergel, Phys. Rev. D 63, 063002 (2001) 

D. Babich and M. Zaldarriaga, Phys. Rev. D 70, 083005 (2004) [arXiv:astro-ph/0408455l . 

M. Liguori, A. Yadav, F. K. Hansen, E. Komatsu, S. Matarrese and B. Wandelt, Phys. Rev. D 76, 
105016 (2007) [Erratum-ibid. D 77, 029902 (2008)] |arXiv: 0708 .37861 [astro-ph]]. 

E. Sefusatti, M. Liguori, A. P. S. Yadav, M. G. Jackson and E. Pajer, JCAP 0912, 022 (2009) 
[arXiv:0906.0232l [astro-ph.CO]]. 

D. Baumann et al. [CMBPol Study Team Collaboration], AIP Conf. Proc. 1141, 10 (2009) 
[arXiv:0811.39l9l [astro-ph]]. 



14 



[29] N. Aghanim, S. Majumdar and J. Silk, Rept. Prog. Phys. 71, 066902 (2008) |arXiv:0711.05l8l [astro- 
ph]]. 

D. M. Goldberg and D. N. Spergel, Phys. Rev. D 59, 103002 (1999) | jarXiv:astro-ph/981125l| . 

D. Hanson, K. M. Smith, A. Challinor and M. Liguori, Phys. Rev. D 80, 083004 (2009) 
[arXiv:0905.4732l [astro-ph.CO]]. 

A. Mangilli and L. Verde, Phys. Rev. D 80, 123007 (2009) [arXiv: 0906 .23171 [astro-ph.CO]]. 
A. Lewis, A. Challinor and D. Hanson, JCAP 1103, 018 (2011) [arXiv:1101.2234l [astro-ph.CO]]. 
P. Serra and A. Cooray, Phys. Rev. D 77, 107305 (2008) |arXiv:0801.3276l [astro-ph]]. 



K. M. Smith and M. Zaldarriaga, arXiv:astro-ph/0612571 



N. Bartolo, S. Matarrese and A. Riotto, Advances in Astronomy, 157079 (2010) 
doi:10.1155/2010/157079, arXiv: 1001.3957vl. 

P. Creminelli and M. Zaldarriaga, Phys. Rev. D 70, 083532 (2004) [arXiv:astro-ph/0405428l . 



N. Bartolo, S. Matarrese and A. Riotto, JCAP 0606, 024 (2006) [arXiv:astro-ph/0604416| . 

N. Bartolo, S. Matarrese and A. Riotto, JCAP 0701, 019 (2007) |arXiv:astro-ph/0610110~l . 

C. Pitrou, Class. Quant. Grav. 24, 6127 (2007) [arXiv:0706.4383 [gr-qc]]; C. Pitrou, Class. Quant. 
Grav. 26, 065006 (20091. larXiv:0809.3036l [gr-qc]. 

L. Senatore, S. Tassev and M. Zaldarriaga, JCAP 0908, 031 (2009) jarXiv:0812. 36521 [astro-ph]]. 
M. Beneke and C. Fidler, Phys. Rev. D 82 (2010) 063509 [arXiv: 1003. 18341 [astro-ph.CO]]. 



C. Pitrou, J. P. Uzan and F. Bernardeau. larXiv:0807.0341l [astro-ph]. 
N. Bartolo and A. Riotto, JCAP 0903, 017 (2009) |arXiv:0811.4584l [astro-ph]]. 

D. Nitta, E. Komatsu, N. Bartolo, S. Matarrese and A. Riotto, JCAP 0905, 014 (2009) 
[arXiv:0903.0894l [astro-ph.CO]]. 

C. Pitrou, J. P. Uzan and F. Bernardeau, JCAP 1007, 003 (2010) |arXiv:1003.048T1 [astro-ph.CO]]. 

E. Komatsu, Class. Quant. Grav. 27, 124010 (2010) |arXiv:1003.6097l [astro-ph.CO]] . 
Talks at [http : / /benasque. or g / 20 1 1 Astrophysics / 1 

N. Bartolo, S. Matarrese and A. Riotto, JCAP 0508, 010 (2005) ||arXiv:astro-ph/0506410"] . 

L. Boubekeur, P. Creminelli, G. D'Amico, J. Norena and F. Vernizzi, JCAP 0908, 029 (2009) 
[arXiv:0906.0980l [astro-ph.CO]]. 

[51] S. Dodelson, Modern cosmology, Amsterdam, Netherlands: Academic Pr. (2003) 440 pp. 



15 



[52] W. Hu, Phys. Rev. D 62, 043007 (2000) [arXiv:astro-ph/0001303] . 

[53] W. Hu and M. J. White, Astrophys. J. 479, 568 (1997) [ arXiv:aatro-ph/9609079| ■ 

[54] C. Cheung, A. L. Fitzpatrick, J. Kaplan and L. Senatore, JCAP 0802, 021 (2008) [arXiv:0709.0295l 
[hep-th]]. 

[55] N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. Lett. 93, 231301 (2004) [arXiv:astro-ph/0407505] . 



16 



